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ABSTRACT 

We consider the interaction between a relativistic fireball and material assumed to be still 
located just outside the progenitor star. Only a small fraction of the expected mass is sufficient 
to efficiently decelerate the fireball, leading to dissipation of most of its kinetic energy. Since 
the scattering optical depths are still large at distances comparable to the progenitor radius, the 
dissipated energy is trapped in the system, accelerating it to relativistic velocities. The process 
resembles the birth of another fireball at radii R ^ 10^^ cm, not far from the transparency 
radius, and with a starting bulk Lorentz factors Fc ^ 10. As seen in the observer frame, this 
"re-generated" fireball appears collimated within an angle dj — l/Fc. If the central engine 
works intermittently, the funnel can, at least partially, refill and the process can repeat itself. 
We discuss how this idea can help solving some open issues of the more conventional internal 
shock scenario for interpreting the Gamma-Ray Burst properties. 

Key words: gamma rays: bursts. X-rays: general, radiation mechanisms: general. 



1 INTRODUCTION 

The internal/external shock scenario (see e.g. Piran 2004; Meszaros 
2006) is currently the leading model to explain the complex phe- 
nomenology of Gamma-Ray Bursts (GRBs) prompt and afterglow 
emission. Despite it can account for many observed characteristics, 
there are a few open issues and difficulties that this model cannot 
solve, or can accommodate only with some important modifica- 
tions. 

Here we recall some problems of the standard scenario and 
mention some ideas already put forward to account for them. 

• Efficiency I: high T-contrast — In internal shocks only the 
relative kinetic energy of the two colliding shells can be dissi- 
pated. Thus "dynamical" efficiencies of only a few per cent can 
be achieved for colliding shells whose Lorentz factors T differ by 
a factor of order unity. Such efficiency has to include energy dis- 
sipated into randomizing protons, amplifying (or even generating) 
magnetic fields and accelerating emitting leptons. As the emitted 
radiation is produced only by the latter component, it corresponds 
to just a fraction of the dynamical efficiency. This problem, pointed 
out, among others, by Kumar (1999) can be solved by postulating 
contrasts in T much exceeding 100 (Belobodorov 2000; Kobayashi 
& Sari 2001). In these cases the typical Lorentz factor of GRBs 
should thus largely exceed the "canonical" value ~ 100. In this 
case it is difficult to understand how the value of the peak energy 
of the prompt spectrum does not wildly change. 

• Efficiency II: Afterglow/Prompt power ratio — A related in- 
consistency concerns the observed ratios of the bolometric fluence 
originated in the afterglow to that in the prompt phase. Since exter- 



nal shocks are dynamically more efficient than internal ones, such 
ratio is expected to exceed one, contrary to current estimates. The 
problem has been exacerbated by the recent observations by the 
Swift satellite, showing that the X-ray afterglow light curve seen 
after a few hours - thought to be smoothly connected with the end 
of the prompt - comprises a steep early phase. As a consequence, 
the total afterglow energy is less than what postulated before. Will- 
ingale et al. (2007), parameterizing the behavior of the Swift GRB 
X-ray light curves, derived an average X-ray afterglow-to-prompt 
fluence ratio around 10 per cent (see also Zhang et al. 2007). Fur- 
thermore, the very same origin of the early X-ray radiation as pro- 
duced by external shocks is questioned, since its behavior is differ- 
ent from the optical one (e.g. Panaitescu et al. 2006). If the X-ray 
emission does not originate in the afterglow phase, this further re- 
duces the above ratio. 

• The spectral energy correlations — Correlations have been 
found between i) the energy where most of the prompt power is 
emitted (iJpeak) and the isotropic prompt bolometric energetics 
E^,iso (Amati et al. 2002; Amati 2006), and ii) -Epeak and the col- 
limation corrected energetics E^. The slope of the former corre- 

1/2 

lation is -Epcak oc -E^ while the slope of the latter depends on 
the radial profile of the circumburst density. For a homogeneous 
density medium Epcuk oc E'^''' (Ghirlanda et al. 2004), while for 
a wind-like profile in density (oc r~^) the correlation is linear: 
Bpcak oc E^ (Nava et al. 2006; Ghirlanda et al. 2007a). If lin- 
ear, the relation is Lorentz invariant and indicates that different 
GBRs roughly emit the same number of photons at the peak (i.e. 
E^/Epcak ~ constant). 

Note that the derivation of E^ requires not only information on 
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the jet break time fj, but also a model relating with the colli- 
mation angle 9j , which in turn depends on the circumburst density 
value, profile and the radiative efficiency -q (i.e. E~i — r;i5kin). The 
phenomenological connection among the three observables E-y^iso, 
i?pcak and tj, as found by Liang & Zhang (2005), is instead model- 
independent. It is of the form EiBo oc E^^^^f^^, which for 6 ~ — 1 
is consistent with the Ghirlanda relation (both in the homogeneous 
and wind case; see Nava et al. 2006). A further tight phenomeno- 
logical relation appears to link three prompt emission quantities: 
the isotropic peak luminosity Liso, Speak and the time interval 
To. 45 during which the emission is above a certain level (Firmani 
et al. 2006). All these correlations were not predicted by the inter- 
nal/external shock scenario, and can only be reconciled with it as 
long as specific dependences of the bulk Lorentz factor upon Eiso 
are satisfied (see Table 1 in Zhang & Meszaros 2002). 

The above issues motivate the search for alternatives or 
for substantial modifications of the standard model. The effi- 
ciency problem and the existence of the spectral-energy rela- 
tions prompted Thompson (2006, T06 hereafter) and Thompson 
Meszaros & Rees (2007, T07 hereafter) to suggest that, besides 
internal shocks, dissipation might also occur because of the inter- 
action between the fireball and the walls of the funnel in the star 
through which it propagates (see §2). This hypothesis also intro- 
duces a typical scale to the problem, namely the radius of the pro- 
genitor star (_R* ^ lO^^-lO^^ cm): shear instabilities within _R, 
can reconvert a significant fraction of bulk kinetic energy into heat. 
Since this dissipation occurs up to _R, (i.e. not far from the trans- 
parency radius), the increased internal energy can only partially re- 
convert into bulk motion via adiabatic expansion, increasing the ef- 
ficiency. Similarly, studies of magnetized fireballs (e.g. Drenkhahn 
& Spruit, 2002; Giannios & Spruit 2007) have shown that dissipa- 
tion of magnetic energy through reconnection can also contribute 
to increase the radiation content of the fireball at relatively large 
radii. 

Along the above mentioned lines, in this Letter we propose a 
further possible way in which a large fraction of the fireball bulk 
energy can be dissipated at distances i? ~ i?, . This assumes that 
at these distances the fireball collides with some mass which is 
(nearly) at rest: a small fraction of the mass swept up in the funnel 
left along the fireball propagation axis is sufficient to lead to effi- 
cient dissipation. Hereafter such mass will be referred to as "IDM" 
(Intervening Debris of the cocoon Material). 

The treatment of the collision is simplified, in order to allow an 
analytical and simple description. We assume the IDM to be at rest 
and homogeneous in density and the fireball to have a Lorentz fac- 
tor r 1. The interaction is described in "steps", while in reality 
will be continuous in time. A complete treatment of the dynamics 
and emission properties of our model requires numerical simula- 
tions (of the kind presented in Morsony et al. 2007, introducing 
some erratic behavior of the injected jet energy). Interestingly, the 
model predicts that jet properties depend on the polar angle (like 
in a structured jet, see Rossi et al. 2002) and it naturally implies a 
connection among the observed spectral-energy correlations. 



2 SHEAR-DRIVEN INSTABILITIES AND DISSIPATION 
OF BULK KINETIC ENERGY 

As mentioned, T06 and T07 proposed a model in which the ef- 
ficiency of the dissipation of kinetic energy into radiation is en- 
hanced with respect to the internal shock scenario and the spectral 
energy correlations, in particular the Amati one, can be accounted 



for. At the same time, in their scenario synchrotron emission could 
play a minor role, the radiation field being dominated by thermal- 
ized high energy photons or by the inverse Compton process. 

For what follows it is useful to summarize their main argu- 
ments here. Consider a fireball that at some distance Rq ~ i?* from 
the central engine is moving relativistically with a bulk Lorentz fac- 
tor Fq. The fireball is initially propagating inside the funnel of the 
progenitor star. T06 and T07 assume that a large fraction of the en- 
ergy dissipated at Ro is thermalized into black-body radiation of 
luminosity 



iBB.iso = 4nRlroaTo''' = 4:-K—^aTQ, 

i n 



(1) 



where Tq and To = FoTq are the temperatures at Ro in the comov- 
ing and observing frame, respectively. The coUimation corrected 
luminosity is Lbb = (1 — cos6'j)LBB,iso which, for small semi- 
aperture angles 9j of the jetted fireball (assumed conical), gives 

2 2I/BB 



(2) 



A key assumption of the model is that Fq ~ 1 /6j . The argument 
behind it is that if Fo > > 1/Oj, shear-driven instabilities have not 
time to grow (in the comoving frame), while in the opposite case 
the flow mixes easily with the heavier material and decelerates to 
F ~ 1. Then, assuming Fo — l/Oj and substituting it in Eq.[T] 



r o d2 ^BB rr 4 
i^BB.iso ~ onUo- alo . 



(3) 



Setting Ebbabo = -^^BB.isotburst and Ebb = ^BBiburst, where 
tburst is the duration of the prompt emission, gives 



rr ,-1-1/4,-1/4 
Speak OC To CX -Bb'b^j^^Sbb' ij^^^^. 



(4) 



This corresponds to the Amati relation if Ebb is similar in differ- 
ent bursts and the dispersion of GRB durations is also limited. It 
should be noted that a relation similar to the Amati one can be also 
recovered adopting -Ebb oc -Ep^ak' suggested by the Ghirlanda 
relation. For instance, for a = 1 (wind case) 



-Epeak OC E-Jg -^^^ty^^^^^. 



-1/5 



(5) 



For the derivation of Eq. 4 (and 5) a key assumption is the depen- 
dence on temperature of the black-body law, which leads to both a 
slope and a normalization similar to those characterizing the Amati 
relation. 

We can ask what happens if, instead of a black-body, one as- 
sumes that the spectrum is a cutoff power-law. This question is par- 
ticularly relevant since the burst spectrum is rarely described by a 
pure black-body (even if some bursts is, see Ghirlanda et al. 2003), 
and also the black-body plus power-law model (Ryde 2005) faces 
severe problems, even considering time resolved spectra (Ghirlanda 
et al. 2007b). 

How is the above derivation modified if, instead of a black- 
body emission, the spectrum is best described by a cutoff power- 
law? Consider then a spectrum described in the comoving frame by 
L'-fAsaiE') OC £■ ^'^ exp(— S'/Sf',) and approximate the observed 
isotropic bolometric luminosity as 



Tiso Fo 



2 ( Speak \ 



1-/3 



-(l+/3)/2 



This leads to 

Speak OC 



2/3).-l/{l-/3) 



^pe'ai (6) 



(V) 
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Therefore the dependence of the Amati relation can be recovered 
even for cutoff power-law spectra, but the normalization in this 
case is not determined. Note also that for [3 — —3 (Wien spectrum), 
Eq.|7]and Eq.|4]have the same dependences. 



3 FIREBALL-IDM COLLISION 

To excavate a funnel inside a progenitor star of mass Af « — 10M« , i 
solar masses, the "proto-jet" has to push out the mass that did not 
fall into the newly born black hole. This means a fraction of (1 — 
cos6»f)M. = 0.1A/.,i6l|_„i, where Oi = Cl^/,-! is the funnel 
opening angle. This mass expands sideways as the proto-jet breaks 
out at the surface of the progenitor, forming a cocoon (see also 
Ramirez-Ruiz, Celotti & Rees 2002). We must expect, however, 
that after the break-out the region in front of the funnel will not 
be perfectly cleared of mass. To be irrelevant, the mass Mc left as 
IDM should be < Eo/{roc^) ~ 5 x lO^^A'/©: less than one in a 
million particles should remain there. 

If the jetted fireball is not continuous, this mass may be still 
there at the moment of the arrival of the new fireball pulse. And 
even if Mc is (10~"^-10~*) the excavated mass, the IDM can have 
important dynamical effects on it. As the bulk velocity and energy 
content of the IDM can be neglected in comparison to those of the 
coming fireball, it will be approximated as initially at rest and cold. 

An interesting aspect of this scenario concerns multi-peaked 
bursts, especially when pulses in the prompt emission are separated 
by quiescent periods. If the central engine works at a reduced rate 
during quiescence, material from the walls of the funnel, previously 
in pressure equilibrium with the jet, will tend to refill again the 
funnel (see also Wang & Meszaros 2007) . This requires a time 
Tf ~ e{R/{l3sc) ~ 3 X 10-*(i?/i?o,7)6lf,_i//3s,-i s, where /3sC is 
the sound speed, and Ro = 10^i?o,7 cm is the radius at the base of 
the funnel. Then the amount of mass which should be pushed out 
again depends upon the quiescent time, but that after ~1 second, 
the funnel is completely closed, and the process repeats itself. 



3.1 Results: dynamics and dissipation 

In the following we derive the main characteristics of the sys- 
tem formed by a fireball impacting against the IDM. Consider a 
fireball with energy Eq, mass A/q and bulk Lorentz factor Fq = 
Eo/ {Mqc?), impacting against a mass Afc, initially at rest. 
The energy and momentum conservation laws read 



100 



MoVo + Mc = rc{Mo + Mc + € /c^) 
MoToPo = r^Pc {Mo + Mc + e'/c'' ) 



(8) 



where e' is the dynamical efficiency measured in the frame moving 
at /3cC. Solving for /3c and e' gives 



/3c 



(3o 



AfoTo 



MoTo + Mc 



Eo 

Tc 

Fee'. 



1 + x ~ xFc 



/3o 
To 



McC^ 
MoToc^ 



(9) 

(10) 
(11) 



The Lorentz factor Fc is shown in Fig. [T] as a function of A/c for 
four values of Eq (top panel) and as a function of Eo for three 
values of Mc (bottom panel). One can see that for the process to 
be interesting (i.e. Fc significantly smaller than Fo) and to avoid 
"over-loading" of baryons (too small Fc) Mc is required to be in 
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Figure 1. The Lorentz factor Fc of the IDM+fireball system as a function 
of Mc for different values of Eq (top panel) and as a function of for 
selected values of Mc (bottom panel). 



specific ranges, that depend on Eo - These ranges, however, encom- 
pass almost two orders of magnitude. 

Fc as a function of x is reported in Fig. |2] The power-law 
dependence Fc oc x~^^'^ can be derived directly from Eq.|9] since 
F7^ = 1 - /3c oc 2x for x < 1 and (io 1- Tc is limited to 
~ 10 even for a:: ~ 4 x 10"'', corresponding to Mc = 2.2 x 

lO~^So,52A/0. 

In Fig.|2]we show t]' = e' /Eo and rj = e/Eo SLsa function of 
X. For clarity also the fraction of the initial kinetic energy preserved 
after the collision 



??kii 



l)(Afo + Afc)c' 
Eo 



= (Fc 



(12) 



is plotted. The sum rj + 77kin is unity by definition. 

Note that since a; is a ratio, Afc and Eo can be taken both as 
"isotropic" or "real" (collimation corrected) values. 



3.2 Evolution of the flreball-i-IDM system 

After the collision/dissipation phase the fireball-l-IDM is expected 
to be optically thick: the Thompson scattering optical depth of the 
IDM material is 



— acuAR — 



OtMc 

AnRinir: 



= 3.2 X 10'' 



Mc. 



-"•c.ll 



(13) 



while the optical depth of the fireball just before the collision is of 
order 



ro = 



(JtEo 



47ri?2ror 



1.8 X 10' 



E-p' 



Rii 



(14) 



where Mc and Eo are here isotropic quantities. These large optical 
depths imply that the radiation produced following the collision is 
trapped inside the fireball-l-IDM system, which will expand because 
of the internal pressure. In the frame moving with Fc, the expansion 
is isotropic. In this frame some final F' will be reached. As seen 
in the observer frame, the expansion is highly asymmetric, and the 
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Figure 2. The Lorentz factor rc,2 = Fc/lOO, the dynamical efficiency e' 
and e (in the comoving and observer frame, respectively) as a function of a;. 
Each quantity is calculated for Fq = 50, 100 and 200 (from top to bottom). 



geometry of the system resembles a cone, with semi-aperture angle 
given by (Barbiellini, Celotti & Longo 2003) 

/3± P'smO' ^ 1 

^ ~ rT' 



T^, (15) 



' rc(/3'cose'+/3c) 

where the last equality assumed 6' = 90° (/3' ~ 1 and /3c ~ 
1). Therefore the aperture angle of the re-born fireball is related to 
Fc, independently of the initial F-factor of the fireball before the 
collision or after the expansion. Note that the initial aperture angle 
of the fireball is irrelevant, as it is the aperture angle of the funnel 
of the progenitor star. 

The fireball is not coUimated in a perfect cone, and mass and 
energy propagate also outside 9j. Since in the frame moving with 
Fc it expands isotropically, M' {yi') = M' / (47r) is approximately 
constant. Therefore in the observer frame 
M' dcosO' 



M{e) 



(16) 



47r d cos 6 

Also the resulting Lorentz factor F is angle dependent: from the rel 
ativistic composition of velocity (e.g. Rybicki & Lightman 1979): 

P' cos e' + /3c 



Pi- 
rn 



/3sin6l = 



(1 + /3c/3' cos 6') 

13' sine' 
Fc(l + /3c/3' cos 61') 

,2 \l/2 



(17) 



+/31, 

F(e) = [l-/3'(e)]~'/' = F'Fc(l + /3'/3cCose'), 

The observed F-factor is constant up to angles slightly smaller than 
l/Fc, and decreases as above. 

As a consequence of the angular dependence of mass and bulk 
Lorentz factor, also the energy depends on 9 as 



E{e) = r{e)M{e)c^ 



(18) 



Such dependences of mass, F and E on polar angle are illustrated 
in Fig. [3] The jet is structured and well approximated by a top-hat 



Figure 3. The mass, energy and bulk Lorentz factor (observer frame) as 
functions of the angle from the jet axis 6, for Fc = 10 and F' = 10. 
All the three quantities are constant up to ~ 1/Fc and then decrease 
(approximately) as power-laws. 



jet: the energy profile is nearly constant within an angle slightly 
smaller than 1/Fc, and at larger angles decreases approximately as 
a steep power-law E{e) oc 6*""/^], since M{e) oc 6"'^^^ and 
r{d) oc for 9 S> l/Fc. This particular behavior gives raise to 
an afterglow light curve indistinguishable from a top-hat jet (see 
Rossi et al. 2004). 



4 DISCUSSION AND CONCLUSIONS 

The most appealing feature of the proposed model is the high ef- 
ficiency in re-converting the fireball kinetic energy into internal 
pressure at a radius comparable to the radius of the progenitor star, 
i.e. on a scale not far from the transparency one. Also the ob- 
served energetics of the internal radiation will be large, since the 
system becomes transparent during (or slightly after) the expan- 
sion/acceleration phase, similarly to the standard fireball models of 
initially high entropy, where the fireball becomes transparent be- 
fore coasting. Our model therefore increases the parameter space 
of high efficiency regimes. Our model is also similar to the model 
proposed by T06 and T07 and to those in which the dissipation 
of an energetically important magnetic field occurs at large radii 
(see e.g. Giannios & Spruit 2007). The efficiency of the energy re- 
conversion for the fireball-IDM collision is of the order of 50-80 
per cent for large ranges in the mass of the IDM and energy of the 
fireball (see Fig.[2ll. 

With respect to the idea proposed by T07, summarised in §2, 
in our scenario there is no requirement on any specific value for 
the fireball bulk Lorentz factor Fo prior to its collision with the 
IDM. Note also that in the T07 model, a "standard" fireball (i.e. not 
magnetic) moving with Fo ~ 1 /9j can dissipate part of its kinetic 
energy, but it cannot reach final F-factors larger than Fo, which is 
bound to be small for typical Oj . In our case the final F can instead 
be large, (even if always smaller than Fo). The jet angle is ~ l/Fc 
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even for large values of the initial To and final F. This is a result of 
our model, and not an assumption. 

The re-born fireball is structured and the M{9) and T{6) be- 
haviors imply that E{9) depends on 9 as steep power-laws. De- 
spite the angle dependence of the energy, the jet should produce an 
afterglow indistinguishable from a top-hat jet. Clearly our descrip- 
tion of the fireball-IDM interaction is extremely simplified, aimed 
at building a physical intuition based on the analytical treatment. 
More realistic situations should be studied via numerical simula- 
tions, but we would like to comment on two aspects, i) Even if the 
IDM is initially at rest, as soon as the fireball starts depositing a 
fraction of energy and momentum, the IDM would begin to move. 
The whole process will take long enough that towards the end it 
would be probably better described as the interaction with a mov- 
ing IDM, with a consequent loss of efficiency, ii) In a 'continuos' 
(non-intermittent) scenario, the IDM would predominantly inter- 
act with the fireball edge, causing only a partial dissipation of its 
energy and the formation of a "fast spine-slow layer" structure. In 
this case the determination of the relevant jet opening angle (i.e. 
within which most of the energy is concentrated) requires a more 
accurate numerical treatment. 

Despite of these caveats, it is still interesting to consider 
whether the model can account for the spectral-energy correla- 
tions, or at least highlight the relations between them. Consider the 
Ghirlanda correlation in the wind case, iJpcak oc E^, that can be 
rewritten as -Epoak oc OfEj^iao for small 9j. The requirement that 

1/2 

also -Bpeak OC E^\^^ (Amatl relation) leads to: 

E-y9f = constant, (19) 

i.e. more energetic bursts are more coUimated, as predicted in our 
model. The above condition (Eq.|19t can be quantitatively satisfied 
if i) the mass against which the fireball collides is similar in dif- 
ferent GRBs, namely Fc ~ l/Sj ~ (see Fig.[B and ii) the 
prompt emission luminosity E^ is also a constant fraction of Eq for 
different GRBs. Eq.[T9]does not explain the Amati or Ghirlanda re- 
lations, although it offers some physical meanings of required con- 
nection between the two. Note that, in the standard internal shock 
model model, one can recover the Amati relation if F ~ constant 
(see Zhang & Meszaros 2002). 

In this Letter we did not discuss the characteristics of the 
spectrum predicted in our scenario (Nava et al., in prep). In gen- 
eral terms, the most effective radiation process would be "dynam- 
ical Compton", or Fermi "acceleration" of photons, as discussed 
by Gruzinov & Meszaros (2000) in the context of internal shocks. 
The large number of photons per proton in the fireball (correspond- 
ing to the "fossil" radiation that has accelerated the fireball itself 
in the first place) implies that a significant fraction of the internal 
energy following the fireball-IDM collision directly energizes pho- 
tons, i.e. photons amplify their energy by interacting with leptons 
with bulk momentum not yet randomized in the shock. The process 
is analogous to particle acceleration in shocks and gives raise to 
high energy photons, conserving their number. 

Since this occurs at i? ~ 10^^ cm, i.e. slightly above the 
progenitor star, the re-bom fireball will become transparent dur- 
ing or just after the acceleration phase (at a transparency radius 
Rt ~ 3 X 10^^ cm, Eq.|I3l). This ensures that photons do not have 
time to loose energy via adiabatic expansion, once again leading to 
large efficiencies. 

The timescale for the refilling of the funnel during quiescent 
phases of the central engine can be short enough that the process 
repeats itself. The estimated refilling timescale is of the order of 



a second: bursts with shorter or longer 'quiescent' phases should 
then show different properties. If the fireball-IDM collision is the 
dominant process for the dissipation (and emission), more ener- 
getic spikes are expected to follow longer quiescent phases (see 
Ramirez-Ruiz & Merloni, 2001). 
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